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CURVE FITS OF PREDICTED INVISCBD STAGNATION- POINT RADIATIVE 


HEATING RATES, COOLING FACTORS, AND SHOCK STANDOFF 
DISTANCES FOR HYPERBOLIC EARTH ENTRY 

By John T. Suttles, Edward M. Sullivan, 
and Stephen B. Margolis 
Langley Research Center 

SUMMARY 

Curve-fit formulas are presented for the stagnation- point radiative heating rate, 
cooling factor, and shock standoff distance for inviscid flow over blunt bodies at condi- 
tions corresponding to high-speed earth entry. The data which were curve fitted were 
calculated by using a technique which utilizes a one-strip integral method and a detailed 
nongray radiation model to generate a radiatively coupled flow-field solution for air in 
chemical and local thermodynamic equilibrium. The range of free-stream parameters 
considered were altitudes from about 55 to 70 km and velocities from about 11 to 
16 km/sec. Spherical bodies with nose radii from 30 to 450 cm and elliptical bodies 
with major -to- minor axis ratios of 2, 4, and 6 were treated. 

Power-law formulas are proposed and a least-squares logarithmic fit is used to 
evaluate the constants. It is shown that the data can be described in this manner with an 
average deviation of about 3 percent (or less) and a maximum deviation of about 10 per- 
cent (or less). The curve-fit formulas provide an effective and economic means for 
making preliminary design studies for situations involving high-speed earth entry. 

INTRODUCTION 

One of the problems encountered by hypervelocity heat-shield designers and trajec- 
tory analysts is that of calculating the radiative flux to an entry body. For economic rea- 
sons, it is desirable to conduct studies and preliminary design work by using simplified 
approaches and correlation equations such as those in reference 1. The radiative heating- 
rate correlation used in reference 1 was derived from early heating-rate predictions 
which were based on a transparent, constant property shock layer. More recent studies 
have shown the early work to be grossly in error for high-speed entry problems. The 
survey given in reference 2 shows that for high-speed entries a heating- rate analysis 
must include the combined effects of shock- layer radiative cooling (which results in 


properties varying across the shock layer), nongray self-absorption, continuum radiation, 
and atomic line radiation to be valid. 

Numerous analyses are now available (e.g., refs. 3 to 8) which include these impor- 
tant effects for inviscid, stagnation region flow. Several of these techniques (refs. 4 to 6) 
have used approximate absorption-coefficient models for the radiation calculations in the 
interest of making engineering simplifications to this complex problem. However, these 
approaches are not in a convenient form for designers and analysts who desire estimates 
of shock- layer radiative heat transfer without resorting to lengthy computer solutions. 
Based on the success of the graphical correlation in reference 9, it is reasonable to 
expect that a simple analytic expression for radiative heating rates, such as that used 
in reference 1, can be derived by curve fitting the results of detailed computer solutions. 
Therefore, such an effort was undertaken. 

In addition to heating rates an effort was made to curve fit the cooling factor since 
it is believed that the results presented in this form will be of more lasting value than the 
heating-rate results. The cooling factor is simply the ratio of the heating rate for a non- 
adiabatic shock layer to the heating rate for an adiabatic shock layer with the heating rates 
based on the same radiation model, body size, and flight conditions. (It is well known that 
for hypervelocity flows nonadiabatic effects significantly reduce the radiative heating com- 
pared to the heating calculated by assuming adiabatic conditions.) 

Olstad (refs. 4 and 9) has demonstrated that the cooling factor can be correlated as 
a function of the heating rate for an adiabatic shock layer normalized by the free- stream 
kinetic- energy flux ^ p The twofold usefulness of the cooling factor was also sug- 

gested by Olstad (ref. 9). First, the use of the cooling factor minimizes the effects of the 
radiation model used in the calculation; heating rates which have significantly different 
values as a consequence of the use of different radiation models show good agreement 
when compared on the cooling-factor basis. This rationale is used for the comparison 
of results in the present report. Second, because the cooling factor has this property of 
bringing radiative heating rates calculated with diverse radiation models to a common 
base, it can be used in conjunction with an adiabatic heating rate calculated with the best 
available radiation model to calculate the nonadiabatic heating rate. With this approach 
it is not necessary to use time-consuming coupled flow-field solutions to obtain radiative 
heating rates, nor is it necessary to recalculate the cooling factors as improvements are 
made in radiation models. 

Since inviscid shock- standoff -distance values are obtained when heating rates and 
cooling factors are calculated, a curve fit for this quantity was derived also. The shock 
standoff distance is of interest to heat-shield designers, for example, in evaluating the 
effects of a probe extending from the body into the shock layer. Also, flow-field analysts 


find that the shock standoff distance has numerous uses, two of which are for comparing 
various results and for initiating iterative solutions. 

The present report presents the results of curve fitting information on calculated 
heating rate, cooling factor, and shock standoff distance with simple expressions involving 
free-stream conditions and vehicle geometry parameters. The calculations were made 
by using a radiatively coupled flow-field computer program (ref. 3) which includes a 
detailed nongray radiation model (refs. 10 and 11). The range of free-stream parameters 
considered were altitudes from about 55 to 70 km and velocities from about 11 to 
16 km/sec. Spherical bodies with nose radii from 30 to 450 cm and elliptical bodies 
with major-to-minor axis ratios of 2, 4, and 6 were treated. The simple form used for 
the curve fits is suitable for hand calculation with a slide rule or electronic desk calcu- 
lator and can be readily included in a parametric analysis. 

SYMBOLS 

ajj coefficient matrix of equation (5) 

Aj constants used in heating-rate curve fit (see eq. (1)) 

bjj coefficient matrix of equation (8) 

Bi constants used in cooling-factor curve fit (see eq. (7)) 

c^ coefficient matrix of equation (10) ' 

constants used in shock-standoff-distance curve fit (see eq. (9)) 

Dj constant vector of equation (5) 

e base for natural logarithms (e = 2.71828 . . .) 

E total-error measure (see eq. (4)) 


F c cooling factor, ratio of nonadiabatic-to-adiabatic heating rates 

Gj constant vector of equation (8) 

Hj constant vector of equation (10) 
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i index used for constants in curve-fit equations 

j index used for equations resulting from minimization procedure 

K constant in functional form assumed for curve-fit equations 

n index used to identify individual data points which are to be curve fitted 

(case number) 

N total number of data points for a curve fit 

q radiative heating rate, W/cm2 

R 

r axis ratio of elliptical bodies, free -stream velocity always oriented normal 

to major axis 

Rg body nose radius, cm 

free-stream velocity, km/sec 
X denotes any of quantities to be curve fitted 

6 stagnation-point shock standoff distance, cm 

e n error of curve fit at the nth data point 

p free-stream density (obtained from 1962 standard atmosphere for a given 

OO 

altitude), g/cm^ 

Subscripts: 

C curve fit 

D data 

r axis ratio 

6 shock standoff distance 
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Superscripts: 


a,b,c exponents in the functional form assumed for curve -fit equations 

METHOD 

The data to be curve fitted were obtained from solutions of the inviscid radiating 
flow field in the stagnation region of a blunt body (fig. 1), which is entering the Earth’s 
atmosphere. These solutions were obtained by using the one- strip integral method of 
reference 3 which incorporates the nongray radiation model of reference 10 as a sub- 
routine. The radiation model ("RATRAP”) is described in further detail and compared 
with other models in reference 11. The computational method is for the inviscid flow of 
air in chemical and local thermodynamic equilibrium. The radiation calculations assume 
a plane slab shock-layer geometry and neglect absorption by the air upstream of the shock 
layer and byproducts of ablation at the body surface. Discussion of the importance of 
these effects can be found in reference 2. 

The data on which the curve fits are based are given in table I. The free- stream 
conditions, velocity and density (density determined from 1962 standard atmosphere for 
a given altitude), were picked to cover the range of interest for high-speed entry cases 
where radiative heating will be important. Figure 2 illustrates the conditions chosen: 
a typical manned planetary -return lifting-entry trajectory (ref. 12), a range of high- 
velocity ballistic -entry trajectories (from an unpublished study), and a matrix of cases 
which span the trajectories. The body sizes were selected from results of an unpublished 
study and cover a range of interest from small unmanned probes of preliminary flight- 
test vehicles to large manned spacecraft. Thus, the curve fits should be of value for a 
wide range of applications. 

The effort was undertaken on the premise that if sufficient data were available it 
would be possible to write a curve -fit equation of the form 

X = Kp a R T D b V c 

^oo ID oo 

where X represents either the heating rate, cooling factor, or the standoff distance and 
a, b, and c may be allowed to vary with the free-stream conditions and/or body size. 
This form has been used to correlate convective heating rates (see, e.g., ref. 13), and it 
is well suited for calculations made with either a slide rule or electronic desk calculator. 
It can also be readily programed for repetitive calculations on a digital computer. 
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CURVE -FIT EQUATIONS 


Heating Rates 

Spherical bodies.- A curve fit of the stagnation- point radiative heating rate q 

rt 

with the free-stream velocity V^, free-stream density p , and the body radius R R 

was obtained by fitting a multidimensional curve through the spherical body data (cases 1 

to 91) in table I. Based on a preliminary analysis, the functional form 


*R,C 


= e 


A 2 +A 3 V oor b A 4 +A 5 V co+A 6 V 00 2^ 


A 7 +A 8 V oo +A 9 VJ 


( 1 ) 


was assumed and the constants Aj, A 2 , . . Ag were found by a least-squares analy- 
sis. To facilitate application of the least-squares method, equation (1) was first linear- 
ized with respect to the A^ values by taking logarithms as follows: 


ln q R,C = A 1 + A 2 ln P oo + A 3 V °° ln p oo + A 4 ln R B + A 5 V oo ln R B 

+ A r V 2 ln R n + A 7 ln + A ft V ln + A Q V 2 ln (2) 

6 °o B 1 00 o °o 00 g 00 00 ' ' 

The error function for the least-squares method is, therefore, 

e n = In q R D (n) - ln q R c (n) (3) 

where q R g(n) is the value of q R corresponding to the nth data point jT^n), V co (n), 

and Rg(n)] from table I, and q R ^(n) is the q R value obtained from the curve fit 
(eq. (1)) for the nth point. The total-error measure in the method is given by 

N 

E - £ *„ 2 ' < 4 > 
n=l 

and is to be minimized with respect to each of the A^ values. The minimization leads 
to a set of linear simultaneous equations for the Aj values whose matrix form is as 
follows: 


a ji A i = Dj (i,j = 1,2,. . ., 9) (5) 

where the A^ values have previously been defined and the coefficient matrix a^ and 
Dj values are given in table n. The elements of the coefficient matrix involve summa- 
tions of functions involving p (n), R R (n), and V (n); and the D. values involve summa 

00 W J 
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tions of functions involving P oo (n), Rg(n) , V^n), and ^(n). In each case, the sum- 
mations are over the N data points. Equation (5) was solved by inverting the coefficient 
matrix a^. The results for the values are given in the appendix, and the results 
using equation (1) with the values are given in table in(a) for the 91 points used. A 
comparison of q^ g with q^ ^ in table in(a) shows an average deviation of approxi- 
mately 3 percent. ’Also, of the 91 data points only 3 have deviations larger than 7 percent, 
and these are for small bodies at high altitudes - conditions for which the calculated 
radiative heat- transfer rates tend to be very low and the inaccuracy may be large. 

Elliptical bodies .- The curve fits presented were based on the data for spherical 
bodies given in table I. Data are given in table I (cases 102 to 150) for elliptical bodies 
with axis ratios of 2, 4, and 6. For the elliptical shapes the free -stream velocity is 
normal to the semimajor axis and Rg is the stagnation- point radius of curvature. 

When the curve fit for the spherical-body data was applied to the nonspherical-body data, 
a systematic error appeared which depended on the axis ratio. This error was attributed 
to the fact that the axis ratio influences the tangential velocity gradient and the shock 
standoff distance. It was found possible to compensate for this effect by use of a correc- 
tion factor (called the shape factor) which is a linear function of the axis ratio. The shape 
factor was obtained by calculating q^ ^ for each of the nonspherical-body cases by 
using equation (1) and forming the ratio of q R ^ to q R g. These ratios were collected 
into groups of data corresponding to common axis ratios, and simple averages for each 

C 

group were determined. The value of = — for the spherical bodies (unit axis ratio) 

4 R,D 

was taken to be 1.0. Thus, values for the shape factor were generated for axis ratios of 
1, 2, 4, and 6. A linear function was assumed for the variation of the shape factor with 
axis ratio, and a least-squares curve-fit method was used to find the slope and intercept. 
Figure 3 presents the data used in the shape-factor calculations and the line obtained 
from the least-squares fit. With the resulting shape factor, the expression for q^ q 
for nonspherical bodies is 


q R,C 

q R,C,r " (0.955859 + 0.03645r) 


( 6 ) 


which is used for 1 < r S 6. It is noted that the shape factor is not precisely 1.0 at unit 
axis ratio since this condition was used merely as a data point rather than as a constraint 
for the curve fit. 

The q,-, values are compared with the corresponding q D n values in 

KjL ,r S\y±J 

table m(b). The results are excellent with an average deviation of 3 percent and a max- 
imum of 8 percent for the cases where Rg, p , and are in the range for which the 

A ^ values were originally calculated. 
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Cooling Factors 


A curve fit for the cooling factor F^. was derived from the 53 cases in table I for 
which Fq values were given. The curve-fit formula for F^ as a function of p , 

V , and Rt 3 was obtained in the same manner as that for q„. A function of the form 

v P B 1_ b 2 +B 3 R B +B 4 R B 2+B 5 R B 3 r B 6 v B 7 n) 

F C,C = e p oo R B V oo (7) 

was assumed based on a preliminary data analysis. Taking the logarithm of equation (7) 
and applying the least- squares method resulted in a set of equations for the B^ values 
whose matrix form is as follows: 

bj^i = G. (ij = 1,2,..., 7) (8) 

where the B^ values are the constants appearing in equation (7) and the coefficient 
matrix b^ and values are given in table IV. The elements of the coefficient 
matrix involve summations over the N data points of functions involving p^(n), V^n), 
and Rg(n), whereas the Gj values involve summations of functions involving p w ( n), 
V^n), Rg(n), and Fq g( n )- Equation (8) was solved by inverting the coefficient matrix 
bjj. The results for the Bj values are presented in the appendix, and the results using 
equation (7) with the Bj values to compute F^-. q are given in table V. A comparison 
of Fq q with F^ g is also given. These results for the cooling factor have an aver- 
age deviation of approximately 3 percent and a maximum deviation of about 10 percent. 

Shock Standoff Distances 

Spherical bodies .- A curve -fit formula for the shock standoff distance at the stag- 
nation streamline as a function of p , V , and R n was obtained in the same manner 

oo 00 D 

as those for q^ and F^. Ninety-one spherical-body cases (cases 1 to 91) were used. 

A function of the form 

C< C 9 +C,V^ C.+CcV^ C fi +C 7 V +C ft V 2 
6 C = e Ip 2 3 Rg 4 5 6 7 « 8 °° ( 9 ) 

was assumed based on a preliminary analysis of the available standoff-distance data. 
Again, by taking the logarithm of the assumed function and applying the least- squares 
method, a set of equations for the values was obtained. The matrix form of these 
equations is as follows: 

c jjCj = Hj (i,j = 1 , 2 ,..., 8) (10) 
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where the Cj values are the constants appearing in equation (9) and the coefficient 
matrix Cj^ and Hj values are given in table VI. The elements of the coefficient 
matrix involve summations over the N data points of functions involving p^(n) , Rg(n), 
and V oc (n); whereas the Hj values involve summations of functions involving P oo (n), 
Rg(n), V TO (n), and 6g(n). Equation (10) was solved by inverting the matrix Cjp and the 
values are given in the appendix. Table VII(a) gives a comparison of the original 
standoff- distance data 6g with 6^ values calculated by using equation (9) and the 
derived values. The results show that when compared to 6g the corresponding 
©C values have an average deviation of about 1 percent with a maximum deviation of 
about 5 percent. 

Elliptical bodies .- A shape factor for the standoff distance for the nonspherical- 
body cases was derived in the same manner as that for the heating- rate correlation. The 
resultant expression for the shock standoff distance for the elliptical bodies is 

6 ?C 

C ’ r (0.825784 + 0.14816r) 

which is applied for 1 < r i 6. The data from which the shape factor for 
was derived are shown in figure 4. The comparison of 6 q r with 6g r 
table VII(b) shows that the curve-fit values have an average deviation of 4 
pared to the original data. 


(ID 

equation (11) 
given in 
percent com- 


DISCUSSION 

As shown in tables El, V, and VE the curve-fit equations give good representations 
of the basic data (table I) in the ranges specified in the appendix. In this section the 
basic data and/or the curve fits will be compared with other available data. In addition, 
the cooling-factor curve fit was studied extensively in order to establish its sensitivity to 
out-of-range conditions. The results of this study are also presented. 

Heating Rates 

Radiative heat-transfer calculations are known to be significantly influenced by both 
the flow-field calculation technique and the radiation transfer model. The basic data 
(table I) for heating rate and standoff distance were calculated by using the one-strip 
integral technique of reference 3, which includes the radiation model of reference 10. 
Results from calculations using these techniques have been compared with several other 
methods and these comparisons have been published (e.g., refs. 3, 7, and 8). The general 
finding is that the one- strip integral technique produces results which compare very well 
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with other techniques, provided that the same radiation model is used. Thus, it can be 
concluded that the fluid mechanics of the method are adequate. 

The effects resulting from the use of different radiation models have also been 
explored (refs. 3 and 11). It has been found that the calculated heating rates for a given 
set of conditions can differ by 50 percent or more. This result is not surprising when 
one considers the diversity in the amount of detail programed into various radiation 
models. In view of the lack of experimental data, at the conditions of interest here, it is 
difficult to assess which of the available radiation models gives the most accurate results 
and which provides the best compromise between accuracy and required computer time. 
For these reasons the present heating-rate results are not compared with previously 
published values. Instead, such comparisons will be made on a cooling-factor basis in 
the next section. 


Cooling Factors 

As indicated in the "Introduction," Olstad (ref. 9) has shown that when the cooling 
factor Fq is taken as a basis, good agreement is obtained between results calculated 
by using different radiation models. Therefore, the F^ curve fit of the present report 
will be compared to the charts of F^ presented by Olstad in reference 9. It is empha- 
sized that to derive heating rates from Fq values a companion calculation of the adia- 
batic radiative heating rate is required. That is, 


( q R) . " F c( q R) 

' 'nonadiabatic ' Adiabatic 


( 12 ) 


and /q R \ does not require a time-consuming coupled flow-field solution. Since 

' ^'adiabatic 

Fq is relatively insensitive to the radiation model (ref, 9), this allows the user to select 
the model most appropriate to his purpose to determine the adiabatic value. If the user 
desires to accept the radiation model applied in the calculation of the basic data of this 
report, he can simply use the heating- rate curve fit presented herein. 


The basic data have been curve fitted wherever possible (the 53 cases having Fq 
values in table I) in terms of Fq; the resulting curve-fit equation is given in the appen- 
dix. It can be seen (table V) that the agreement between the basic data and the curve fit 
is excellent with a maximum error of approximately 10 percent. As an independent check, 
a large number (180) of cooling-factor values were calculated by using the curve-fit equa- 
tion (eq. (7)) and compared with the curves given in reference 9. Figures 5, 6, and 7 show 
the results of these calculations. Note that some calculations shown in the figures are 
out of the range of the correlation. This was done to show the sensitivity of the curve fit. 
(The calculations also included cases with Rg = 800 cm; however, they produced an 
unreasonably small Fq and have not been plotted.) In figure 5 it is seen that the values 
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from the curve fit for those cases where all three parameters are within the ranges spec- 
ified in the appendix are in excellent agreement with the results of reference 9. The 
maximum difference is approximately 17 percent for one case of a 30-cm body at an 
altitude near the upper limit of the altitude range. The values from figures 5, 6, and 7 
for cases where one or more of the parameters are out of range show that the curve fit 
is most sensitive to body size and least sensitive to velocity. 

In view of the results shown in figures 5, 6, and 7, the authors believe that the 
cooling-factor curve fit is accurate provided the ranges specified in the appendix are not 
exceeded. 


Shock Standoff Distances 

ft 

Figure 8 presents a comparison of the shock standoff distances from the basic data 
(table I), the curve-fit equation (the appendix), and the results of Callis (ref. 6) for 
p^ = 1.225 x 10-7 g/cm3 and for = 15.24 km/sec and 12.19 km/sec. The basic 
data do not contain identically corresponding results, and, therefore, data for the closest 
conditions ^p = 1.2959 x 10"7 g/cm3 and = 15 km/sec and 12 km/secj are given. 

It is seen from figure 8 that these slight variations in and p^ have little effect 

on the comparison. The comparison shows that the curve-fit equation gives an excellent 
representation of the basic data and, although the trends are different, the basic data and 
the results of Callis agree within 10 percent or less. A more detailed comparison 
between the two methods (i.e . , Suttles (ref. 3) and Callis (ref. 6)) is given in reference 3, 
where it is shown that the differences between radiation models can account for a signifi- 
cant difference in the results. Therefore, the variation in the trends of the results in 
figure 8 are most probably due to the different radiation models used. 

Some comparisons were made (but not shown here) between the results of the pres- 
ent work and the results of Callis for bodies smaller ^Rg = 1 to 30 cmj and for bodies 
larger ^Rg = 450 to 1000 cmj than those considered in the curve fits. This comparison 
indicated that for bodies in the range of Rg from 10 to 30 cm, the results remained in 
good agreement (10 percent or less); but for smaller bodies and bodies larger than 450 cm, 
the results began to deviate significantly. It is concluded that, as in the case of the 
cooling-factor curve fit, the user should be very hesitant in applying the shock- standoff- 
distance curve fit outside the range of the basic data. 

, CONCLUDING REMARKS 

Curve-fit formulas are presented for the stagnation- point radiative heating rate, 
cooling factor, and shock standoff distance for inviscid flow over blunt bodies at condi- 
tions corresponding to high-speed earth entry. The data on which the curve fits are 
based were calculated by using a technique which utilizes a one-strip integral method 
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and a detailed nongray radiation model to generate a radiatively coupled flow-field solu- 
tion for air in chemical and local thermodynamic equilibrium. The range of free-stream 
parameters considered were altitudes from about 55 to 70 km and velocities from about 
11 to 16 km/sec. Spherical bodies with nose radii of from 30 to 450 cm and elliptical 
bodies with major-to-minor axis ratios of 2, 4, and 6 were treated. 

Power -law formulas are proposed and a least- squares logarithmic fit is used to 
evaluate the constants. It is shown that the data can be described in this manner with 
an average deviation of about 3 percent (or less) and a maximum deviation of about 
10 percent. A study of the sensitivity of the formulas indicates that they should be used 
only within the range of the free stream and vehicle geometry parameters of the data. 
These curve-fit formulas provide an effective and economic means for making prelimi- 
nary design studies for situations involving high-speed earth entry. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 30, 1974. 
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APPENDIX 


CURVE -FIT FORMULAS AND RELATED PARAMETERS 

For the convenience of the reader, the curve-fit formulas, the constants derived 
for the formulas, the ranges of applicability of the formulas, and the average deviation 
of the formulas from the basic data are listed in this appendix. 

Heating Rates 

Formula .- The formula for in W/cm2 is 

n . _ K A l n A 2 +A 3 V °o R A 4 +A 5 v °°+ A 6 V oo 2 v A 7 +A 8 V 00 +A 9 V 00 2 
q R,C ” K r e P oc R B y oo 

Derived constants. - 


r i 

(0.955859 + 0.03645r) 
. 1 


(r*l) 
(r = 1) 


Aj = -69.099 Ag = 0.005381 

A 2 = 1.320 Arj = 51.89 

Ag = -0.01223 A g = -1.558 

A 4 = 1.688 A 9 = 0.02659 

A 5 = -0.1796 

Ranges of applicability . - 

= 1.078 X 10-7 to 6.53 x 10-7 g/ cm 3 
Altitude = 53.75 to 68.4 km 

Rg = 30 to 450 cm 

\ 

= 11 to 16 km/sec 
r = 1 to 6 x 
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APPENDIX 


Average deviation .- 3 percent 

Cooling Factors 

Formula.- The cooling factor F^ q is dimensionless and is given by 


C,C 


B 1 B 2 +B 3 R B +B 4 R B 2+B 5 R B 3 „ B 6„ B 7 


= e p 


r B V co 


Derived constants. - 


B l = -3.679 


B 5 = 0.000000005567 


B 2 = -0.3598 


Bg = 1.059 
B ? = -1.931 


B 3 = 0.002024 
B 4 = -0.000005583 
Ranges of applicability . - 

= 5.61 x 10-7 to 8.75 x 10'8 g/ cm 3 
Altitude = 55 to 70 km 
Rg = 30 to 450 cm 
= 11 to 16 km/sec 
Average deviation .- 3 percent 

Shock Standoff Distances 

Formula.- The formula for shock standoff distance in cm is 


6 C = K r,6 e P 


Derived constants. - 


C 0 +C,V_ C /1 +C c V„ c 6 +c 7 v..+c q v^ 2 


l n ^2 3 v °o R '~4‘ rv " / 5 v oo 


B 


'8 00 


K r,6 = 


(0.825784 + 0.14816r) 
1 


(r * 1) 
(r = 1) 
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APPENDIX 


C 1 = -3.697 
C 2 = 0.07375 
C 3 = -0.003338 
C 4 = 1.134 

Ranges of applicability . - 
p = 1.078 xl0- 7 1 

OO 


Cg = -0.01248 
C g = 0.7718 
C 7 = -0.02572 
C g = 0.00009347 

6.53 x 10'^ g/cm3 


Altitude = 53.75 to 68.4 km 
Rg = 30 to 450 cm 
V M = 11 to 16 km/sec 
r = 1 to 6 


Average deviation.- 2 percent 
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TABLE I.- BASIC DATA 


Case 

V 

v oo> 

km/sec 

Alt., 

km 

Poo’ 

g/cm3 

R B’ 

cm 

Axis 

ratio 

6 d> 

cm 

q R,D’ 

W/cm2 

f c,d 

i 

11.000 

38.0 

3 . 9072 E -07 

30.00 

1.0 

1.421 

306.0 

- 

2 

11.000 

58.0 

3 . 9072 E -07 

150.00 

1.0 

7 . C 36 

544.0 

- 

3 

11. 000 

58.0 

3 . 9072 fc - 0 7 

300.00 

1.0 

13.974 

66 1.0 

- 


11.000 

58.0 

3 . 90726-07 

450.00 

1.0 

20.844 

752.0 

- 

5 

13. 000 

68.0 

3.90722 07 

30.00 

1.0 

1.36 8 

1607.0 

.582 

o 

13.000 

58.0 

3 . 90726-07 

150.00 

1.0 

6.680 

2263.0 

- 

7 

13. 000 

58.0 

3 . 90726-07 

300.00 

1.0 

12.887 

2699.0 

- 

6 

13.000 

58.0 

3 . 90726-07 

450.00 

1.0 

18.971 

3080.0 

- 

9 

14.000 

58.0 

3 . 9072 E -07 

30.00 

1 .0 

1.339 

2578.0 

.508 

10 

14.000 

58.0 

3 . 90726-07 

150.00 

1.0 

6 . 379 

3537.0 

- 

U 

14.000 

58.0 

3. 90 726-07 

300.00 

1.0 

12. 268 

4289.0 

.307 

u 

1 «*. 000 

58.0 

3.90 72 E -07 

450.00 

1.0 

17.943 

4816.0 

.285 

13 

15. 000 

58.0 

3 • 90 72 E - 07 

30.00 

1.0 

1. 30 8 

3709.0 

.452 

14 

15.000 

58.0 

3 . 9072 E - C 7 

150.00 

1.0 

6 . 093 

5124.0 

.298 

15 

15.000 

58.0 

3 . 9072 E - C 7 

300.00 

1.0 

11.611 

6108.0 

.268 

16 

15.000 

58.0 

3. 90 72 E- C 7 

450.00 

1.0 

16. 850 

6917.0 

.252 

17 

16. 000 

58.0 

3 . 9072 E -07 

30.00 

1.0 

1.277 

5462.0 

.401 

18 

16.000 

58.0 

3.9072 6-07 

150.00 

1.0 

5.830 

7039.0 

.261 

19 

16. 000 

58.0 

3 • 90 72 E - 0 7 

300.00 

1.0 

11.011 

8176.0 

.231 

20 

16.000 

58.0 

3 . 9072 E - C 7 

450 . 0 C 

1.0 

15.873 

9382.0 

.224 

21 

11. 000 

61.0 

2 . 7018 E - 07 

30.00 

1.0 

1. 397 

194.0 

.812 

22 

11.000 

61.0 

2 • 70 18 E - 07 

150.00 

1.0 

6.921 

365.0 

.626 

23 

11.000 

61.0 

2 . 7018 E - 07 

300.00 

1.0 

13.725 

447.0 

.541 

24 

11.000 

61.0 

2 . 70 18 E -07 

450.00 

1.0 

20.464 

509.0 

.509 

25 

12.000 

61.0 

2 . 70 18 t - 07 

30.00 

1.0 

1.378 

520.0 

.710 

26 

12.000 

61.0 

2 . 7018 E -07 

150.00 

1.0 

6 . 760 

843.0 

.497 

27 

12. 000 

61.0 

2 . 7018 E -07 

300.00 

1.0 

13.390 

958.0 

.412 

28 

12.000 

61.0 

2 . 7018 E - C 7 

450.00 

1.0 

19.844 

1087.0 

.383 

29 

13.000 

61.0 

2 . 7018 E -07 

30.00 

1.0 

1.348 

1027.0 

.630 

30 

13.000 

61.0 

2 . 7018 E - 07 

150.00 

1.0 

6.640 

1545.0 

.429 

31 

13.000 

61.0 

2 • 7 Q 18 E -07 

300.00 

1.0 

12.710 

1786.0 

.361 

32 

13.000 

61.0 

2 . 7018 E - 07 

450.00 

1.0 

18.800 

1951.0 

.327 

33 

14.000 

61.0 

2 . 7018 E - C 7 

30.00 

1.0 

1.319 

1737.0 

.566 

34 

14.000 

61.0 

2 . 70 18 E -07 

150.00 

1.0 

6.280 

2353.0 

.356 

35 

14.000 

61.0 

2 . 7018 E -07 

300.00 

1.0 

12.090 

2839.0 

.318 

36 

14.000 

61.0 

2 . 7018 E -07 

450.00 

1.0 

17.660 

3243.0 

.303 

37 

15.000 

61.0 

2 . 70186-07 

30.00 

1.0 

1.287 

2656.0 

.514 

38 

15.000 

61.0 

2 . 7018 E - 07 

150.00 

1.0 

5.995 

3471.0 

.314 

39 

15.000 

61.0 

2 . 7018 E - C 7 

300.00 

1.0 

11.440 

4154.0 

• 282 

40 

15.000 

61.0 

2 . 7018 E - C 7 

450.00 

1.0 

16.850 

4690.0 

.268 
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TABLE I.- BASIC DATA - Continued 


Case 

Voo, 

km/sec 

Alt., 

km 

P oo’ 

g/cm 3 

r b> 

cm 

Axis 

ratio 

<5d> 

cm 

q R,D > 

W/cm2 

f c,d 

41 

16.000 

6 1.0 

2 . 7018 E - C 7 

30.00 

1.0 

1.256 

3778.0 

.466 

42 

16.000 

61.0 

2 . 7018 E -07 

150.00 

1.0 

5.713 

4771.0 

.272 

43 

16.000 

61.0 

2 . 7018 E -07 

300.00 

1.0 

10.920 

5362.0 

.233 

44 

16.000 

61.0 

2 . 7018 E -07 

450.00 

1.0 

15.550 

6174.0 

.227 

45 

11.000 

67.0 

1 . 2959 E - C 7 

30.00 

1.0 

1.360 

70.0 

- 

46 

U .000 

67.0 

1. 2959 E — 07 

150.00 

1.0 

6.754 

149.0 

- 

47 

11.000 

67.0 

1 . 2959 E - C 7 

3 CO . 00 

1.0 

13.445 

190.0 

- 

48 

11.000 

67.0 

1 . 2959 E ~ 07 

450.00 

1.0 

20.019 

215.0 

- 

49 

12.000 

67.0 

1 . 2959 E - C 7 

30.00 

1.0 

1.342 

194.0 

- 

50 

12.000 

67.0 

1 . 2959 £-07 

150.00 

1.0 

6.597 

363.0 

- 

51 

12.000 

67.0 

1 . 2959 E -07 

300.00 

1.0 

13.045 

442.0 

- 

52 

12.000 

67.0 

1 . 2959 E -07 

450.00 

1.0 

19. 392 

490.0 

- 

53 

13.000 

67.0 

1 . 2959 E -07 

30.00 

1.0 

1.317 

39 o .0 

- 

54 

13.000 

67.0 

1 . 2959 E -07 

150.00 

1.0 

6.392 

674.0 

- 

55 

13.000 

67.0 

1 . 2959 E -07 

300.00 

1.0 

12.525 

796.0 

- 

56 

13.000 

67.0 

1 . 2959 E -07 

450.00 

1.0 

18.498 

860.0 

- 

57 

14.000 

67.0 

1 . 2959 E -07 

30.00 

1.0 

1.292 

688.0 

- 

58 

14.000 

67.0 

1 . 2959 E -07 

150.00 

1.0 

6 . 194 

1079,0 

- 

59 

14.000 

6 7.0 

1 . 29596-07 

3 C 0.00 

1.0 

12.022 

1218.0 

- 

60 

14.000 

67.0 

1 . 2959 E -07 

450.00 

1.0 

17.636 

1347.0 

- 

61 

15.000 

67 . 0 

1 . 2959 E - C 7 

30.00 

1.0 

1.268 

1064.0 

.620 

62 

15.000 

67.0 

1 . 2959 E - 07 

150.00 

1.0 

5.980 

1576.0 

.362 

63 

15.000 

67.0 

1.2959 E - 07 

300.00 

1.0 

11.470 

1750.0 

.291 

64 

15.000 

67.0 

1 • 2959 E - C 7 

450.00 

1.0 

16. 710 

1913.0 

.269 

65 

16.000 

67.0 

1 . 2959 E -07 

30.00 

1.0 

1.242 

1531.0 

- 

66 

16. 000 

67.0 

1 . 2959 E-C 7 

150.00 

1.0 

5.756 

2150.0 

- 

67 

16.000 

67.0 

1 • 2959 E - 07 

300.00 

1.0 

10.920 

2360.0 

- 

68 

16.000 

67.0 

1 . 2959 E -07 

450.00 

1.0 

14.850 

2450.0 

- 

69 

13.000 

55.0 

5 . 6080 E - C 7 

30.00 

1.0 

1.394 

2388.0 

- 

70 

13.000 

55.0 

5 . 6080 E -07 ■ 

150.00 

1.0 

6.725 

3411.0 

- 

71 

13.000 

55.0 

5 • 6080 E -07 

300.00 

1.0 

13.064 

4156.0 

- 

72 

13. 000 

55.0 

5 . 6080 E -07 

450.00 

1.0 

19. 194 

4752.0 

- 

73 

15. 000 

55.0 

5 . 6080 E - C 7 

30.00 

1.0 

1.332 

5851.0 

- 

74 

15.000 

55.0 

5. 6080 E -07 

150.00 

1.0 

6.209 

7632.0 

.291 

75 

15.000 

55.0 

5 . 6080 E -07 

3 CO • 00 

1.0 

11.770 

9381.0 

.268 

76 

15.000 

55.0 

5 . 6080 E — 07 

450.00 

1.0 

17.040 

10877.0 

.255 

77 

16. 000 

55.0 

5 . 6080 E - C 7 

30.00 

1.0 

1.298 

8247.0 

- 

78 

16.000 

55.0 

5. 6080 E -07 

150.00 

1.0 

5.895 

10294.0 

- 

79 

16.000 

55.0 

5 . 6080 E - C 7 

300.00 

1.0 

11.070 

12360.0 

- 

80 

16.000 

55.0 

5 . 6080 E — 07 

450.00 

1.0 

15.918 

14314.0 

- 
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TABLE I.- BASIC DATA - Continued 


Case 

V 

v ooj 

km/sec 

Alt., 

km 

P CO’ 

g/cm3 

R B’ 

cm 

Axis 

ratio 

6 D’ 

cm 

q R,D’ 

W/cm2 

f c,d 

81 

15.025 

68. 3 

1.0780E-07 

342.70 

1.0 

12.978 

1542.0 

- 

82 

14.897 

66.3 

1.<V270E-07 . 

342.70 

1.0 

13.040 

2036.0 

- 

83 

14.737 

64. 7 

1.7580E-07 

342.70 

1.0 

13.243 

2348.0 

- 

84 

14. 549 

6 3.4 

2.01 70E-07 

342.70 

1.0 

13. 383 

2616.0 

- 

85 

14.346 

62.9 

2. 1600E-C7 

342.70 

L.O 

13.535 

2618.0 

- 

86 

14.161 

62.8 

2. 1680E-07 

342.70 

1.0 

13.669 

2435.0 

- 

87 

13.560 

62. 8 

2. 1680E-07 

342.70 

1.0 

14.068 

1892.0 

- 

88 

13.650 

63.5 

2. 0259E-07 

60.96 

1.0 

2.634 

1210.0 

- 

. 89 

13. 169 

59. 3 

3.3412E-C7 

60.96 

1.0 

2.716 

1732.0 

- 

90 

12.845 

57.4 

<f.2147E-07 

60.96 

1.0 

2.766 

1883.0 

- 

91 

12.009 

53.7 

6.5255E-C7 

60.96 

1.0 

2.904 

1771.0 

- 

92 

15. 000 

59.5 

3.2535 E-07 

150.00 

1.0 

6. 079 

4187.0 

.304 

93 

15. 000 

59.5 

3.2535E-C7 

3 CO . 00 

1.0 

11.608 

4936.0 

.268 

94 

15.000 

59.5 

3. 2535 E-07 

450.00 

1.0 

16.870 

5531.0 

.252 

95 

15.000 

64.0 

1.8837E-07 

150.00 

1.0 

6.015 

2311.0 

.328 

96 

15.000 

64.0 

1.8837E-07 

300.00 

1.0 

11.530 

2582.0 

.272 

97 

15.000 

64.0 

1. 8837E-07 

450-00 

1.0 

16.625 

2854.0 

.253 

98 

15.000 

70.0 

8.7535E-08 

30.00 

1.0 

1.258 

635.0 

.692 

99 

15.000 

70.0 

8.7535E-08 

150. 00 

1.0 

5.944 

1031.0 

.413 

100 

15.000 

70.0 

8.7535E-08 

300.00 

1.0 

11.468 

1169,0 

.324 

101 

15.000 

70.0 

8. 7535E-C 8 

450.00 

1.0 

16.599 

1252.0 

.289 

102 

13.650 

63.5 

2.02596-07 

121.92 

2.0 

4.688 

1514.0 

- ■ 

103 

13. 169 

59. 3 

3. 3412c- 07 

121.92 

2.0 

4.843 

2090.0 

- 

104 

12.845 

57.4 

4.2147 E-07 

121.92 

2.0 

4.934 

2230.0 

- 

105 

12.458 

55. 5 

5.2997E-C7 

121.92 

2.0 

5.044 

2236.0 

- 

106 

12.009 

53.7 

6. 5255 E- 07 

121.92 

2. 0 

5.173 

2092.0 

- 

107 

11.499 

52.0 

8. 04 13E-0 7 

121.92 

2.0 

5.290 

1723.0 

- 

108 

13.650 

63.5 

2. 02 59E-07 

243.84 

4.0 

7.268 

1639.0 

- 

109 

13.169 

59. 3 

3.3412E-C7 

243.84 

4.0 

7.501 

2210.0 

- 

110 

12.845 

57.4 

4.2147E-07 

243.84 

4.0 

7.625 

2419.0 

- 

111 

12.458 

55.5 

5. 2997E-07 

243.84 

4.0 

7.783 

2480.0 

- 

112 

12.009 

53. 7 

6. 5255E-07 

243.84 

4.0 

7.978 

2325.0 

- 

113 

11.499 

52.0 

8.0413E-07 

243.84 

4.0 

8.242 

1966.0 

- 

114 

13.650 

63.5 

2. 0259 E- 07 

365.76 

6.0 

8.728 

1720.0 

- 

115 

13. 169 

59. 3 

3. 34 12E-07 

365.76 

6.0 

8.963 

2331.0 

- 

116 

12.845 

57.4 

4. 2 1 4 7E-07 

365.76 

6.0 

9. 178 

2553.0 

- 

1 1 7 

12.458 

55.5 

5.2997E-07 

365.76 

6.0 

9. 352 

2628.0 

- 

118 

12. 009 

53. 7 

6. 5255E-07 

365.76 

6.0 

9. 575 

2512.0 

- 

119 

11.499 

52. 0 

8. 0413E-07 

365.76 

6.0 

9.800 

2263.0 

- 

120 

13.580 

64. 2 

l . 844 7E-07 

243.70 

4.0 

6.994 

1375.0 

- 
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TABLE I.- BASIC DATA - Concluded 


Case 

v,, 

km/sec 

Alt., 

km 

Poo’ 

g/cm3 

r b> 

cm 

Axis 

ratio 

6 D’ 

cm 

q R,D’ 

W/cm2 

121 

13.232 

60. 1 

3.0232E-C7 

243.70 

4.0 

7. 195 

2041.0 

122 

12.701 

56.2 

4. 8433E- 07 

243.70 

4.0 

7.423 

2570.0 

1 23 

11.949 

52.6 

7.4351E-C7 

243.70 

4.0 

7. 791 

2565.0 

124 

10.963 

49.3 

1.1199 E— 06 

243.70 

4.0 

8. 134 

1612.0 

125 

12.700 

59.2 

3.3958E-07 

179.83 

4.C 

5.644 

1650.0 

126 

12.254 

55.4 

5.3436E- 07 

179.83 

4.0 

5.817 

2049.0 

127 

11.611 

51.9 

8. 1349 E-07 

179.83 

4.0 

6.063 

1994.0 

128 

13. 161 

61.0 

2.6849E-07 

231.64 

4.0 

7.094 

1800.0 

129 

12.767 

57.1 

4.3343E-C7 

231.64 

4.0 

7.37C 

2365.0 

130 

12. 190 

53.4 

6.7347E-07 

231.64 

4.0 

7.545 

2703.0 

131 

11.398 

50.0 

1. 0247 E— 06 

231.64 

4.0 

7.947 

2339.0 

132 

13.467 

64. 2 

1.8330E-07 

262.12 

4.0 

7.840 

1327.0 

133 

12.909 

59.2 

3.3693 E— 07 

262. 12 

4.0 

8.073 

1964.0 

134 

12.220 

55.5 

5. 2983E-07 

262.12 

4.0 

8.408 

2168.0 

135 

11.295 

52.0 

7.9770E-07 

262.12 

4.0 

8.847 

1635.0 

136 

13.050 

61.1 

2. 666 7 E-07 

252.98 

4.0 

7.724 

1652.0 

137 

12.592 

57.2 

4.2847E-07 

252.98 

4.0 

8.220 

2169.0 

138 

11.932 

53.6 

6.61 35 fc-07 

252.98 

4.0 

8.336 

2255.0 

139 

11.052 

50.3 

9.9419E-07 

252.98 

4.0 

8.695 

1607.0 

140 

12.700 

59.2 

3.3938E-07 

249.94 

4.0 

7.788 

1763.0 

141 

12.254 
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TABLE m.- COMPARISON OF BASIC DATA AND CURVE -FIT RESULTS FOR HEATING RATE 
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TABLE HI.- COMPARISON OF BASIC DATA AND CURVE-FIT RESULTS FOR HEATING RATE - Continued 
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TABLE III.- COMPARISON OF BASIC DATA AND CURVE-FIT RESULTS FOR HEATING RATE - Continued 
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TABLE IV.- MATRIX ELEMENTS AND CONSTANT VECTOR FOR COOLING-FACTOR CURVE FIT 
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rOO'O'OO^OOaOvO^NQOU'NOvOO'OfnO^^^ >0 CT 

5_ C ooaoor^orofnoh->tr-'h*r*>ooo(^^f\4^o>tinr«-p^0'^ 

o ^ no®»ON>nn>OHfNjooopgvONco' 0 ' 0 'r^o'mHa'oo 

^ ||9I| || t | l|(|f|ltf|f|||||l 

U % C\J -* | fONH O 43 O ^<3 rO 

Wg, *i i » 7> iiltitl 

I 


O ooginooNN^aJMn^in^^^^^^pg^m^pnfvr^ 

- cooo^g-o^in®^fnM(roooN«n^<vo^cg>o^OO 

U inmmcg>tnfNiNm(MNf>j®inu>^N>t't>t40^mfninpri 

•••••••••••♦•••••••*•**••• 


n NQ0^vAfNJ00C0rg-4-H^^cg>0^CT*Oh-Nf0Oa'^r»O^ 

- cooooom<^'4jtf>o«orncg*Mf\44-o^cr-^ao^n<S<jrs|'Oin 

O ininro{N^'rgcMrsi^rgc^(NQo^minf^4*^rn>0^“rnrniAm 

h . • . 


* -2 o © 

3 -h • • 

3 2 


OOOOOOOOOOOOOOOOOOOOOOOOOO 


1 M («jl ^ r-4 <-4 



OOOOOOOOO 
ooooooooo 
• •••••••• 

OOOOOOOOO 
rnmoinroir>oin<n 
fO ^ H ff) ^ 


OOOOOOOO 

OOOOOOOO 

• # ♦ * 9 * 9 * 

OOOOOOOO 
inPini^ioomi’i 
h m vf ~4 ra >t 


OOOOOOOOO 

OOOOOOOOO 

OOOOOOOOO 
inomninomnin 
^ rn >t -* cn * 



UOOOOOWOQt^OOQOOOOOOOOOOOOO 
I ( I i I i f i l i I i t i t t I ! I i i i I I I t 

UJUJUJUJJJLUUJUJUJLUUJUJUJUJUJUJUJUJUJllJUiUJUJlU LU LU 

rgfNjogcgogosl^ojrsjrvjrvjfsjaoaooocOaoOQOOaooaoaOcoao 

"■* Hf4-4Hf4w|^HHp4H»4^ 

OOOOOOOOOOOOOOOOOOOOOOOOOO 
O' O' O' O O' O' O' 0 s O'O' ^r"»NN^SNNNNNNf*NS 


rtmrO^rornmrnmrn^rorsjfvjogiNjpgoaosjtNp^rgrgrg^NN 



OOOOOOOOOOOOOOOOOOOOOOOOOO 

aoooaoocDcoaoOoocOcO 

urN\A\AVAtfN\rMnVA^tf\tf%tAs0'0>0^>'0'0N0O'0'0>0 >0 nO NO 



OOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOO 

H H H H r- 4 r*4 »— 4 «*4 —4 H ^4 H -*4 ^ >4 H H *~4 i«4 4 «4 <4 «“l .-4 


a> 

3 inO*-<tNi'n4 , tf>'0^<30CT'O-4fNjr04 , U>>0f^C0l>O-^N^'t 

rt rH^H^4^4-4-4-4-H-4(\trsiCNjogf\|»NrNjfMCMPgmmrO'* f >rri 


28 



TABLE V.- COMPARISON OF BASIC DATA AND CURVE-FIT RESULTS FOR COOLING FACTOR - Concluded 
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TABLE VI.- MATRIX ELEMENTS AND CONSTANT VECTOR FOR STANDOFF-DISTANCE CURVE FIT 
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TABLE VII.- COMPARISON OF BASIC DATA AND CURVE-FIT RESULTS 
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Shock wave 



Figure 1.- Parameters of inviscid radiating flow-field solutions 
for stagnation region of a blunt body. 



Figure 2.- Free-stream conditions and body sizes used in 
generating data to be curve fitted. 
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Axis ratio 

Figure 3.- Shape-factor correction for heating rate of elliptical bodies 
(free stream alined with semiminor axis). 
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Figure 4.- Shape-factor correction for stagnation- point shock standoff distance of 
elliptical bodies (free stream alined with semiminor axis). 
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Curve-fit calculations 

Within range of R g and altitude of present 
curve fit 
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curve fit 

© Altitude too low (48.77 km) for present 
curve fit 


e 



10~ 4 10~ 3 10 ' 2 10' 1 


10 


Adiabatic heating rate / 1 

Figure 5.- Comparison of cooling-factor correlation curves from reference 9 with results 
of present curve fit for three velocities within the range of the present curve fit. 
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velocity = 18.28 km/sec) higher than those of the present curve fit. 




40 


Figure 7.- Comparison of cooling-factor correlation curves from reference 9 with results of present curve fit for 

a velocity = 10.67 km/secj lower than those of the present curve fit. 
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□ Basic data : = 1 5 km/sec, 

feo = 1.2959 X 10 -7 g/cm 3 

- Curve-fit equation: V = 15*24 km/sec, 

feo = 1.225 X 10" 7 g/cm 3 

_ Caliis (ref. 6): = 1 5-2k km/sec, 

ffco = 1 .225 X 10 -7 g/cm 3 
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(b) High velocity. 

Figure 8.- Comparison of shock standoff distances for basic data, 
curve-fit equation, and results of Caliis (ref. 6). 
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